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Microsecond resolution of quasiparticle tunneling in the single-Cooper-pair-transistor 
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We present radio- frequency measurements on a single-Cooper-pair-transistor in which individual 
quasiparticle poisoning events were observed with microsecond temporal resolution. Thermal ac- 
tivation of the quasiparticle dynamics is investigated, and consequently, we are able to determine 
energetics of the poisoning and un-poisoning processes. In particular, we are able to assign an 
effective quasiparticle temperature to parameterize the poisoning rate. 



Operation of both the single-Cooper-pair-transistor 
(SCPT) and the Cooper-pair-box (CPB) rely on the co- 
herent tunneling of a single Cooper-pair between a reser- 
voir and a tunnel-coupled island. This coherent phe- 
nomenon is the basis of CPB charge qubits JSLJ3 and 
low-dissipation electrometry using the SCPT |3,I3- One 
of the challenges that faces these devices is avoiding the 
incoherent tunneling of quasiparticles, often referred to 
as quasiparticle poisoning. The effect of quasiparticle 
poisoning is to change the charge on the device island by 
an electron and halt the coherent tunneling of Cooper- 
pairs. This is especially undesirable for the CPB qubit 
where it can be a major source of decoherence j^^lS- 

Quasiparticle poisoning has been extensively studied 
with a wide variation in behavior observed 4, 7, 8, 9, 10]. 
However, a model suggested by Aumentado et al. ap- 
pears to successfully explain the phenomenon 11]. In 
this model there is some un-known (and possibly non- 
equilibrium) source of quasiparticles in the device leads. 
These quasiparticles are able to tunnel onto (poisoning) 
the device island which usually acts as a quasiparticle 
trap. Subsequently the trapped quasiparticle is ther- 
mally excited (un-poisoning) out of the trap, and the 
island returns to its even-state. 

While most previous investigations of quasiparticle 
poisoning have been performed with a relatively low 
bandwidth, we note a very recent careful study of oxygen- 
doped aluminium SCPTs measured by an rf-technique 
sensitive to the Josephson inductance 12]. In that case, 
detailed measurements of the temperature dependence 
of the poisoned state lifetime allowed determination of 
a quasiparticle trap depth on the island. In this paper 
we present measurements of a SCPT, made by a differ- 
ent technique, embedded in a related radio- frequency (rf) 
tank circuit. A temperature dependent study allowed the 
energetics of both the poisoning and un-poisoning pro- 
cesses to be determined. In particular, a measurement of 
the thermal activation of the un-poisoned state lifetime 
enabled an effective quasiparticle temperature to be de- 
duced which is an experimentally useful way to parame- 
terize quasiparticle poisoning. 

We engineer the SCPT to have a greater supercon- 
ducting gap (2A) for the island than the leads A^ > A^ 
[ll| by making use of the rapid enhancement of A with 
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FIG. 1: (a) Scanning electron micrograph of the device, (b) 
Simplified rf circuit diagram. The rf-carrier signal reflected 
from the tank circuit is amplified by a cryogenic amplifier 
with a gain of 38 dB at 4.2 K. After a further 30 dB of am- 
plification at room temperature, the carrier is phase shifted 
and homodyne detected, (c) Profile of the SCPT showing 
aluminum film thickness and the resulting change in A. 



decreasing film thickness [iS], |l^ (fig. 1(c)). This re- 
duces the depth of the quasiparticle trap, allowing quasi- 
particles to be more easily thermally excited out from 
the island. The island is made from a 10 nm thick film 
(Ai = 250 ±15 /ieV) while the leads have a thickness of 
30 nm (A^ = 208 ± 10 /ieV), with A determined by mea- 
suring the onset of quasiparticle tunneling in SIS junc- 
tions. The device pattern (fig. 1(a)) is defined in polymer 
bilayer resist by electron beam lithography and the alu- 
minum thermally evaporated at a rate of 0.1 nms~^ onto 
a liquid nitrogen cooled stage. With this technique we 
were able to achieve electrically continuous films down 
to a thickness of 5 nm JI3. A controlled oxidization (35 
mTorr for 2 minutes) between the evaporations defines 
the tunnel barriers. 

The circuit used for this experiment is the same as 
commonly used for the rf-SET (fig. 1(b)) jlHl^. A 
resonant circuit at 326 MHz is formed by a chip inductor 
(L = 470 nH), a parasitic capacitance {Cp = 0.51 pF) 
and the SCPT. The reflection coefficient (S'n = \Zt - 
50/Zt+50|) of a small (~ /iV) incident rf-carrier signal at 




FIG. 2: (a) Coulomb diamonds showing the 2e-periodic su- 
percurrent at zero-bias and e-periodic transport at finite bias, 
(b) Amphfied mixer output for Vds = showing the form of 
the supercurrent oscihations. (c) Network analyzer measure- 
ment (heavily averaged) of magnitude and phase of rf-carrier 
across the supercurrent oscillations. The incident power is 
-107 dBm. (d) Model of the amplitude and phase response 
across the supercurrent oscillations. 



the circuit resonance (o;^ = 1/LC) is determined by mis- 
match of the tank circuit impedance {Zt = L/RC) to a 
50 Q co-axial cable. The reflected carrier signal is then 
amplified by a low- noise cryogenic amplifier. Following 
further amplification at room temperature, the rf-carrier 
is demodulated by mixing with a local oscillator at the 
carrier frequency (a technique sensitive to both phase 
and amplitude of the reflected rf-carrier). The resulting 
intermediate frequency (IF) output is further amplified 
and recorded on an oscilloscope. 

When Coulomb diamonds are measured a 2-e peri- 
odic supercurrent is observed at zero bias along with e- 
periodic features at finite bias due to a combination of 
Josephson quasiparticle (JQP) resonances and Coulomb 
blockade of quasiparticle tunneling (fig. 2(a)). Similar 
behavior was seen in a number of other devices. The 
4.2 K resistance of this device was 47 kQ and the charg- 
ing energy Ec = e'^ /2Cy: = 77 /ieV, as determined from 
normal- state Coulomb diamonds measured at 5 = 2.5 T. 
Estimating the Josephson energy per junction from the 
4.2 K resistance and the Ambegoakar-Baratoff relation 
(Ej - 4(A7WS^) ^« fi^d Ej = 33 A^eV. 

Taking a single trace over the supercurrent oscillations 
at Vds = (fig. 2(b)), a change in polarity of the mixer 
output occurs indicating a phase shift of the reflected 
rf-carrier. Further investigation with a network analyzer 
(fig. 2(c)), shows large changes in both the amplitude and 
phase across the supercurrent oscillations. For the am- 
plitude component, there is a high reflection coefficient 
at both rig = 1 (on supercurrent maxima) and n^ = 0, 
and a minima in reflection coefficient on the sides of the 
supercurrent oscillation. There is also a large phase-shift 



{60 = 178 degrees) between rig = 1 and n^ = 0, with the 
phase shifts coinciding with the amplitude minima. 

To understand the behavior of this circuit we develop 
a model in which the rf response depends on the ratio of 
the driving current (J^/ ~ Vrf/ooL = 1 nA at -107 dBm) 
to the switching current (Isw) of the SCPT. We calculate 
switching currents using a 2-band model of the SCPT [3 
finding a maximum Isw =4.6 nA for the ground-band 
at n^ = 1. If Irf < Iswi we assume the SCPT remains 
in current-mode and presents zero resistance. This is 
impedance transformed by the tank circuit to yield Zt ~ 
00 and causes almost complete reflection. 

In the case where Irf > hw^ we assume a hysteresis 
loop at the carrier frequency in which the device is partly 
(for Irf{t) < Isw) in current mode and partly in voltage 
mode (for Irf{t) > Isw)- This leads to an average re- 
sistance {R) which is transformed by the tank circuit to 
Zt = L/{R)C. Using Isw from the 2-band model, an ex- 
pression for (R) 113 and the tank circuit parameters, we 
simulate the amplitude and phase response of the device 
(fig. 2(d)). We note a reduced value of £^j = 8 /ieV was 
taken to account for the suppression of Isw due to envi- 
ronmental effects. The amplitude and phase response are 
well-modeled, with the large phase shift occurring as the 
resonant circuit changes between under {ug = 0, (R) > 1 
Mft) and over-damping {ug = 1, (R) = Q). A phase 
shift is expected in the SCPT due to the Josephson in- 
ductance [131, however in this case it can be explained 
by the resonant circuit going through critical damping. 

Charge sensitivity is determined by applying a 1 MHz 
gate signal of 0.026 e rms and measuring the signal to 
noise ratio (SNR) of the resulting sidebands with a spec- 
trum analyzer. Since there is a phase component to our 
signal we perform this measurement after de-modulation. 
Using the formula Aqrms x 10" "^^^/^^/V^, where B is 
resolution bandwidth, we find a sensitivity of 1.5 x 10~^ 
ellz~^-^ which is comparable to superconducting and 
normal state rf-SETs 17]. 

Two-level switching behavior occurs with the device 
biased on a supercurrent peak and Vif monitored as a 
function of time. The inset in fig. 3(a) shows a represen- 
tative 4 ms time record. The positive level at Vif = 1.6 
V corresponds to the the top of the supercurrent peak 
while the negative level Vif = —1.0 V corresponds to 
the signal in the trough between peaks. We attribute 
the positive voltage state to an even (un-poisoned) state 
where there are only Cooper pairs on the island. By con- 
trast, negative voltages correspond to a 'poisoned-state' 
where a single quasiparticle occupies the island and the 
supercurrent peak is shifted from n^ = 1 to n^ = 0. We 
record traces of length 400 ms, consisting of 10^ data 
points, in order to obtain reliable statistics. 

State parity is ascertained by comparing Vif to a 
threshold half-way between the even and odd levels. The 
distribution of times spent in the even and odd-states is 
measured and plotted in a histogram (fig. 3(a)). We fit 
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FIG. 3: (a) A histogram of times spent in both the even 
and odd-states, the solid lines are fitted exponentials. In- 
set: Switching between even and odd- states observed in non- 
averaged measurements at n^ = 1 with an incident rf-carrier 
power of -108 dBm. (b) The time constants, deduced from 
the previous histogram, as a function of the rf-carrier power. 



an exponential decay e"^/^^ to the histogram, using the 
parameter tj to define the even and odd-state lifetimes. 
For the data in fig. 3(a) todd = 35 /is and teven = 110 
/is. The data is well-fitted indicating that the tunnel 
processes obey Poissonian statistics. A recent study of 
two-level systems has shown that finite receiver band- 
width can have a significant effect on time constants re- 
sulting from an analysis of this type U^. The majority 
of time constants measured are > 10/is and our receiver 
bandwidth is 1 MHz hence, using the analysis in 20], we 
expect the resulting systematic error to be < 10%. 

We measure the even and odd-state lifetimes as a func- 
tion of rf-carrier power to investigate the effect of the 
rf- measurement (fig. 3(b)). A strong reduction in todd 
and only a slight reduction in teven is noticed. The likely 
cause of reduction of the odd-state lifetime is that the rf- 
carrier causes heating which thermally activates a quasi- 
particle off the island. The small reduction in teven indi- 
cates that the rf-carrier doesn't significantly increase the 
quasiparticle population in the leads, and any increase in 
temperature doesn't strongly affect the poisoning rate. 

With the aim of determining the thermal activation 
of the poisoning and un-poisoning events we study the 
even and odd-state lifetimes as a function of tempera- 
ture. A low rf-power (-112 dBm) was chosen to minimize 
heating by the carrier signal. A reduction in todd occurs 
as the temperature is increased, which agrees with the 
measurements in 12j and is due to thermal excitation 
of quasiparticles out of the quasiparticle trap formed on 
the island (fig. 4(a)). Considering the free energy change 
of this transition, we expect the time constant to be ap- 
proximately todd = ^^ ^ Je — ~~ i2ii] 5 where 5E is the 



quasiparticle trap depth, and R the average tunnel junc- 
tion resistance. Fitting the data [23 to this thermal acti- 
vation model we find an experimental value oi5E = 50±4 
/ieV for the trap depth. The expected quasiparticle trap 
depth can be determined by considering the energy differ- 
ence between the poisoned and un-poisoned states |ll|. 
At the supercurrent peak {rig = 1) the maximum un- 
poisoned energy (corresponding to the excited-state, and 
assuming no superconducting phase difference across the 
device) is Eu = Ec -\- ^ ^ while the poisoned state en- 
ergy has a minimum energy of Ep = 5 A = A^ — A/. It 
is useful to note that trap depth changes with gate bias 
(ug) and an investigation of this was carried out in [l^ . 
For our analysis we use the values of Ec = 77 /ieV and 
Ej = 33 /ieV measured for the SCPT and ^A = 42 ± 18 
/ieV as measured from the SIS junctions. Calculating 
SE = Eu — Ep di quasiparticle trap with depth 52 ± 18 
/ieV is expected, hence close agreement to the experi- 
mental value is shown. 

A constant teven is measured up to T ~ 180 mK, 
with the value being reduced by thermal activation at 
higher temperatures. The poisoning rate is expected to 
depend linearly on both the tunnel barrier conductance 
Gi (taken to be in units of ^) and the density of quasi- 
particles in the leads. An expression for the poisoning 
rate was deduced in 6] to be t~J^^^ ^ AT^ exp(— A//A:T) 
for the case where temperature is small compared to trap 

depth (kT « 5E) and A = Qi^^,/^^^. 
The temperature dependence of the poisoning rate can 
be explained by a constant low temperature quasipar- 
ticle population (causing rate t~^) and, at higher tem- 
peratures, the presence of thermally excited quasiparti- 
cles. Adding the rates due to these two populations we 
obtain t'^l^ = t'^ + AT^ exp(-A^/A:T). Figure 4(a) 
shows a fit to the data, with all parameters free, and 
shows good agreement with some deviation above 300 
mK. The parameters we deduce from the fit are a con- 
stant low-temperature even- state lifetime of tc = 102 ± 2 
/is, A = 1.8±0.2 X 10^° K-^-^s-^ and A^ = 213±31 /ieV. 
The value of A^ is in close agreement with the measure- 
ments from the SIS junctions. The calculated value of A 
is 0.92 X 10^^ i^~^"^5~^, hence we see approximate agree- 
ment between the theory and the temperature dependent 
experimental poisoning rate. Poisoning rate (slightly) in- 
creases with trap depth in this regime and the theoretical 
underestimate may be related to the greater trap depth 
indicated by measurement of todd- 

We are able to define an effective quasiparticle tem- 
perature {Tqp). This is the temperature that causes a 
quasiparticle tunneling rate, due to the thermal excita- 
tion of quasiparticles in the leads, equal to the constant 

1 
low-temperature value {t~^ = ATqpexp{— Ai / kTqp)) . For 

this device Tqp = 228 mK, which is slightly greater than 
the electron temperature (Tg ~ 150 mK) as estimated 
from the fit to todd- Due to the relatively long recom- 
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FIG. 4: (a) The time constants teven and todd determined 
as a function of temperature with an rf-carrier power of -112 
dBm. The fit to todd is from a model of thermal activation of 
quasiparticles off the island. For teven, the fitted line includes 
both a constant quasiparticle poisoning rate and thermal ac- 
tivation of quasiparticles across the superconducting gap. (b) 
Even and odd lifetimes as a function of magnetic field. 



bination time of quasiparticles (1 — 10 /is) |23J, quasi- 
particles created by microwave radiation in the leads can 
cause Tqp > T^. For previous measurements (on devices 
with un-engineered A) in which a lack of poisoning was 
observed e.g. 7], either a low quasiparticle temperature 
was achieved or measurement bandwidth was insufficient 
to resolve poisoning events. 

We also perform a quantitative study of the even and 
odd state lifetimes in the presence of an in-plane mag- 
netic field, noting that magnetic fields have previously 
been used to change the periodicity of CPBs [l^ . Little 
change is noticed in teven as the field is increased (fig. 
4(b)). Presumably the poisoning rate remains approxi- 
mately constant until the quasiparticle trap becomes a 
barrier and quasiparticles have to be thermally excited 
onto the island. However, there is a large reduction of 
toddi indicating that the quasiparticle trap becomes shal- 
lower and quasiparticles can more easily be thermally 
excited out. For thinner films the critical field increases 
|l3J . indicating a greater reduction in A^ than A^ (ef- 
fectively increasing 5/S) at a finite fields. For example, 
performing a fit to a temperature dependence of todd at 
B = 150 mT, we find a reduced value of SE = 27 ±2 
/ieV. 

In summary, we employed the change in A with alu- 
minum film thickness to fabricate a SCPT with a reduced 
quasiparticle trap on the island. Individual quasiparticle 
poisoning events were measured and the resulting statis- 
tics analyzed to determine time constants for the even 
and odd-state lifetimes. An important aspect of the ex- 
periment was measurement of the thermal activation of 



poisoning rate. This enabled a quasiparticle tempera- 
ture to be determined which will be a useful parame- 
ter to compare different devices and experimental setups. 
Furthermore, we expect that by reducing the island film 
thickness (increasing ^A) to create a quasiparticle bar- 
rier on the island, or by introducing quasiparticle traps 
(therefore increasing teven) ^ we will be able to fabricate 
SCPTs and CPBs with negligible quasiparticle poisoning. 
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